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Abstract 

Background: The unfolded protein response (UPR) is a network of intracellular signaling pathways that supports 
the ability of the secretory pathway to maintain a balance between the load of proteins entering the endoplasmic 
reticulum (ER) and the protein folding capacity of the ER lumen. Current evidence indicates that several pathogenic 
fungi rely heavily on this pathway for virulence, but there is limited understanding of the mechanisms involved. The 
best known functional output of the UPR is transcriptional upregulation of mRNAs involved in ER homeostasis. 
However, this does not take into account mechanisms of translational regulation that involve differential loading of 
ribosomes onto mRNAs. In this study, a global analysis of transcript-specific translational regulation was performed 
in the pathogenic mold Aspergillus fumigatus to determine the nature and scope of the translational response to ER 
stress. 

Results: ER stress was induced by treating the fungus with dithiothreitol, tunicamycin, or a thermal up-shift. The 
mRNAs were then fractionated on the basis of ribosome occupancy into an under-translated pool (U) and a 
well-translated pool (W). The mRNAs were used to interrogate microarrays and the ratio of the hybridization signal 
(W/U) was used as an indicator of the relative translational efficiency of a mRNA under each condition. The largest 
category of translationally upregulated mRNAs during ER stress encoded proteins involved in translation. 
Components of the ergosterol and GPI anchor biosynthetic pathways also showed increased polysome association, 
suggesting an important role for translational regulation in membrane and cell wall homeostasis. ER stress induced 
limited remodeling of the secretory pathway translatome. However, a select group of transcription factors was 
translationally upregulated, providing a link to subsequent modification of the transcriptome. Finally, we provide 
evidence that one component of the ER stress translatome is a novel mRNA isoform from the yvcl gene that is 
induced by ER stress in a UPR-dependent manner. 

Conclusions: Together, these findings define a core set of mRNAs subject to translational control during the 
adaptive response to acute ER stress in A. fumigatus and reveal a remarkable breadth of functions that are needed 
to resolve ER stress in this organism. 
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Background 

The opportunistic mold pathogen Aspergillus fumigatus 
causes life-threatening pulmonary infections that have 
the potential to progress into invasive aspergillosis, a 
disseminated disease with a very high rate of mortality 
[1,2]. Infections with this fungus continue to impede the 
successful management of patients with hematologic 
malignancies or solid-organ and bone marrow trans- 
plants worldwide, accounting for the highest per person 
hospitalization costs of all the systemic mycoses [3-5]. 
The ongoing expansion of the immunosuppressed po- 
pulation is expected to increase the incidence of the dis- 
ease, which is galvanizing studies to understand more 
about fungal stress response pathways that could yield 
novel vulnerabilities for future therapeutic targeting. 

Current evidence indicates that pathogenic fungi are 
under endoplasmic reticulum (ER) stress in the host 
environment and therefore depend upon adaptive 
stress responses pathways to support their survival 
during infection [6-10]. The unfolded protein response 
(UPR) is the major ER stress response pathway, re- 
sponsible for maintaining an ER lumenal environment 
that is conducive to optimal protein folding [11]. A. 
fumigatus depends upon the UPR to support the expres- 
sion of clinically relevant traits such as thermotolerance, 
cell wall/membrane homeostasis, hypoxia adaptation, iron 
homeostasis, nutrient assimilation from complex sub- 
strates and antifungal drug resistance [6,7]. Similar fin- 
dings have also been reported in Cryptococcus neoformans 
[8], Candida albicans [12], Candida glabrata [10], and 
Alternaria brassicicola [9], suggesting that the UPR is used 
by diverse fungal pathogens as a regulatory hub for the 
expression of multiple attributes that promote virulence 
in the host. The UPR is triggered in response to the accu- 
mulation of unfolded proteins, a condition that arises du- 
ring infection when there is an imbalance between the 
level of nascent proteins entering the ER and the ability of 
the organelle to process that load. ER protein folding may 
also be perturbed by adverse conditions encountered in 
the host such as mammalian body temperature, oxidative 
stress, hypoxia and nutrient limitation [13]. The UPR 
counters the resulting ER stress by expanding the quantity 
of ER-resident chaperones and folding enzymes that are 
needed to help membrane and secreted proteins achieve 
their native conformation. The current understanding of 
the fungal UPR is based upon the paradigm established in 
the model yeast Saccharomyces cerevisiae [14]. The path- 
way is controlled by Irel (IreA in A. fumigatus), an ER- 
transmembrane protein that detects disturbances in the 
ER that lead to the accumulation of unfolded proteins. 
Irel contains a lumenal sensing domain and a cytosolic ef- 
fector region that contains dual enzymes: a kinase linked 
to an endoribonuclease (RNase). In the absence of ER 
stress, Irelp exists as an inactive monomer in association 



with the ER-resident chaperone, Bip (also known as 
Kar2p). When the folding capacity of the ER is ex- 
ceeded, BiP dissociates from the lumenal domain to as- 
sist with protein folding. This triggers the activation of 
Irel by a ligand-dependent two-step mechanism in 
which BiP dissociation is followed by direct interaction 
of Irel with unfolded proteins [15-18]. These events 
elicit Irel oligomerization in the ER membrane, resul- 
ting in a conformational change that activates the 
C-terminal RNase [19,20]. The substrate of this RNAse 
is a cytoplasmic mRNA known as HAC1 (hacA in A. 
fumigatus). The excision of an unconventional intron 
from the HAC1 mRNA allows in-frame translation of 
the bZIP transcription factor, Hacl (HacA in A. fumiga- 
tus). Hacl re-establishes ER homeostasis by remodeling 
the transcriptome to enhance the protein folding capa- 
city of the ER. 

Genome-wide expression profiling has demonstrated 
that A. fumigatus responds to acute ER stress by up- 
regulating the levels of a core group of mRNAs that 
encode proteins with functions that support the se- 
cretory pathway [6]. However, mRNA abundance mea- 
surements do not take translational efficiency into 
consideration, which is a mechanism of gene regula- 
tion that can have potent effects on protein production 
[21-23]. Translational regulation provides the cell with 
a rapid-response mechanism to fine-tune protein levels 
in proportion to need, and is particularly important in 
situations where an immediate response to an environ- 
mental stress is key for survival [24,25]. Translational 
regulation can be studied on a global scale by interro- 
gating microarrays with mRNAs that have been frac- 
tionated based upon ribosome occupancy [26]. This 
approach is based on the fact that translationally quies- 
cent mRNAs are sequestered within messenger ribonu- 
cleoprotein (mRNP) particles or associated with single 
ribosomes (monosomes), whereas actively translated 
mRNAs are associated with multiple ribosomes (poly- 
somes). The hybridization of a microarray with these 
polysome-fractionated mRNAs can thus provide in- 
sight into how the translational efficiency of individual 
mRNAs is modified by environmental cues. Analogous 
approaches have been used to study the ER stress 
translatome in S. cerevisiae and Aspergillus niger 
[27,28]. However, a global analysis of transcript- 
specific translational regulation has not been per- 
formed in A. fumigatus. In this study, polysome 
fractionation of mRNA was coupled with microarray 
detection in order to identify changes in the transla- 
tional status of the A. fumigatus transcriptome under 
conditions that perturb ER homeostasis. The findings 
establish a core ER stress translatome and uncover evi- 
dence for extensive translational regulation during the 
response of A. fumigatus to ER stress. 
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Results & discussion 

Translatome remodeling is a major component of the 
ER stress response in A. fumigatus 

The two most commonly used compounds to induce ER 
stress are dithiothreitol (DTT), which reduces disul- 
fide bonds, and tunicamycin (TM), which inhibits the 
N-linked glycosylation that is required for optimal folding 
[29]. We have previously shown that treatment of A. fumi- 
gatus with 1 mM DTT or 10 ug/ml TM for 1 h is suffi- 
cient to trigger the UPR, defined by the induction of hacA 
splicing and a subsequent remodeling of the transcrip- 
tome to strengthen the protein folding capacity of the ER 
[6,7]. In the present study, we found that these conditions 
had undetectable effects on the overall polysome distribu- 
tion in A. fumigatus (data not shown); both control and 
treated cultures showed a typical polysome profile com- 
prised of individual ribosome subunits, the 80S mono- 
some peak, and polysome peaks representing 2-10 
ribosomes per mRNA (Figure 1). This is similar to S. cere- 
visiae, where treatment with DTT did not result in 
substantial reductions in global translation initiation 
efficiency [27]. However, it contrasts the situation in 
metazoans, where the analogous treatments induce global 
translation attenuation due to phosphorylation of eukar- 
yotic translation initiation factor 2a (eIF2a), thereby redu- 
cing the total burden on the secretory pathway [30]. The 
apparent absence of a global translation attenuation re- 
sponse to ER stress in fungi is consistent with current evi- 
dence that the fungal kingdom lacks the ER stress sensor 
that controls the activation of this pathway [11,31]. 

We hypothesized that a subset of mRNAs that are im- 
portant to surviving ER stress would redistribute into the 
polysome peak to enhance their rate of translation. To test 
this, a genome-wide perspective of mRNA translational 
efficiency during ER stress was obtained by interrogating 
microarrays with mRNAs that were fractionated on the 
basis of ribosome occupancy. ER stress was induced by 
treating the fungus with DTT or TM, as detailed in 
Methods. Ribosome-associated mRNAs were then frac- 
tionated from cytoplasmic extracts into two pools: an 
under-translated pool (fraction-U) containing mRNAs 
with four or less ribosomes and a well-translated pool 
(fraction- W) containing mRNAs with five or more ribo- 
somes (Figure 1). Each fraction was then used to interro- 
gate high-density microarrays and an estimate of the 
translational efficiency of each mRNA was defined here as 
the ratio of the hybridization signal in fraction-W over 
that of fraction-U. Those mRNAs that showed a two-fold 
change (up or down) in this translational efficiency ratio 
during ER stress were considered to be subject to trans- 
lational regulation during ER stress (see Methods for 
additional detail). To maximize the detection of ER stress- 
responsive mRNAs, and minimize chemical-specific ef- 
fects, the resulting dataset was restricted to mRNAs that 



Fraction U 



(under-translated) 
1 




Figure 1 Strategy for investigating the translational efficiency 
of mRNAs by polysome profiling and microarray hybridization. 

A representative polysome profile shows the monosome peak (1) at 
the top of the gradient (left), followed by peaks representing 2,3,4 
and > 5 ribosomes per mRNA. Following centrifugation, the gradient 
was divided into two: an under-translated fraction containing mRNAs 
associated with 1-4 ribosomes (fraction-U) and a well translated fraction 
containing mRNAs associated with 5 or more ribosomes (fraction-W). 
The mRNAs in each fraction were then used to interrogate microarrays, 
as detailed in Methods. The translational efficiency of each mRNA was 
defined as the ratio of the hybridization signal (fraction-W/fraction-U). 
Those mRNAs with a higher W/U ratio during ER stress than in the 
absence of ER stress (using a 2-fold change between conditions as the 
cut-off) were considered to be subject to translational upregulation 
during ER stress. 



showed differential polysome association in response to 
both DTT and TM treatment. Of the 323 mRNAs that fit 
these criteria (Figure 2) the majority showed an increase 
in translational efficiency (233), suggesting that ER stress- 
induced translational regulation is predominandy an in- 
ductive response. It is noteworthy that the number of 
translationally regulated mRNAs in this analysis was over 
four times the size of the transcriptional response pre- 
viously identified under the same conditions [6], demon- 
strating that remodeling of the translatome is a major 
component of the ER stress response in A. fumigatus. 
Validation of the microarray data was obtained by qPCR 
for 3 upregulated mRNAs (ergl, vpsS5, and jlbA) and one 
negative control ((^-tubulin) (Additional file 1). 

A total of 90 mRNAs showed a decrease in transla- 
tional efficiency in response to DTT and TM treatment 
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Figure 2 Translationally regulated transcripts during ER stress 
in A. fumigatus. The Venn diagram demonstrates the number of 
transcripts with a differential translational efficiency ratio in the 
presence of DTT or TM relative to untreated cultures (minimum 2-fold 
change in the W/U ratio). The region of overlap illustrates the number 
of mRNAs that are shared between the two treatments, representing 
mRNAs that are enriched for functions in ER stress response. The 
overlap contains translationally upregulated mRNAs (shaded region) as 
well as downregulated (Down). The area of each circle is scaled to the 
number of transcripts and the values for DTT and TM represent the 
total number of genes that were differentially regulated. 



(Figure 2), 46 of which were unannotated. Analysis of 
the remaining 44 mRNAs using the Kyoto Encyclopedia 
of Genes and Genomes (KEGG) database revealed en- 
richment in metabolism and RNA processing categories, 
consistent with a shift towards more limited metabolic 
functions under ER stress conditions (Additional file 2). 
Of the 1253 mRNAs that showed altered translational 
efficiency in response to DTT alone, 536 were down- 
regulated. This contrasts a study of translational effi- 
ciency during DTT-induced stress in A. niger, which 
revealed down-regulation of 253 mRNAs and upregula- 
tion of 26 mRNAs. A direct comparison between these 
two datasets is difficult however, because the A. niger 
paper treated the fungus with 20 mM DTT for 2 h and 
fractionated mRNAs based upon occupancy with < 2 
and >2 ribosomes [28], whereas our study used 1 mM 
DTT for 1 h and fractionation into pools associated 
with <5 and >5 ribosomes. The more severe ER stress 
conditions used in the A. niger study may account for 
the predominance of translationally repressed mRNAs in 
that organism, relative to the largely inductive response 
in A. fumigatus. 

ER stress induces limited remodeling of the secretory 
pathway translatome 

The transcriptional response of A. fumigatus to acute ER 
stress is narrowly focused on upregulating the level of 
mRNAs that encode proteins that support the secretory 
pathway at multiple levels, including functions such 
as folding, glycosylation, ER-associated degradation, ER 
translocation, vesicular transport and membrane function 



[6]. The majority of these proteins (70%) contained pre- 
dicted signal peptides that would direct them into the ER. 
However, only 18% of the translationally upregulated 
mRNAs identified in the present study contained signal 
peptides. Furthermore, a KEGG pathway- enrichment 
analysis revealed that only a minority (16%) of the mRNAs 
in the translationally upregulated dataset encode proteins 
that directly influence the secretory pathway (Table 1). In 
fact, none of the translationally upregulated mRNAs in 
this study showed increased abundance in our previous 
analysis of the UPR transcriptome under the identical 
conditions [6]. This demonstrates that the nature and 
scope of ER stress-induced translational regulation is fun- 
damentally different from the transcriptional response to 
the same stress. Nevertheless, the identification of several 
transcription factor-encoding mRNAs in the transla- 
tionally upregulated dataset (Table 1) suggests a mechan- 
ism whereby rapid changes in the translatome could be 
linked to subsequent modification of the transcriptome 
during ER stress. Together, these mechanisms would pro- 
vide the fungus with extensive options to modify the 
proteome and restore ER homeostasis. 

ER stress increases the translational state of the 
translation machinery 

We found that components of the translational machinery 
were subject to increased polysome association in the 
presence of either DTT or TM (Table 1). This was some- 
what surprising, since previous studies have shown down- 
regulation of ribosome biogenesis genes in A. niger and 
S. cerevisiae exposed to DTT [27,28]. This discrepancy is 
likely to reflect the higher concentrations of DTT used in 
those studies, and/or species-specific differences in sensi- 
tivity to DTT. We speculate that a limited expansion of 
the translational apparatus is beneficial to A. fumigatus 
during ER stress because it provides a mechanism to 
rapidly increase the level of proteins that are needed to 
protect the ER from damage until the appropriate tran- 
scriptional modifications can be implemented. Since only 
a subset of the translational machinery was upregulated in 
A. fumigatus, a second possibility is that some of these 
proteins may have unrecognized 'moonlighting' functions 
that are relevant to ER stress responses, a possibility that 
is supported by an emerging literature on extra-ribosome 
functions for ribosomal proteins [32] . 

ER stress induces remodeling of the cell wall and 
membrane translatome 

The major interface between A. fumigatus and the host 
environment is the plasma membrane and cell wall, both 
of which are important targets for current antifungal 
therapy [33]. Damage to either of these structures re- 
quires the delivery of new cell wall and membrane com- 
ponents to the hyphal tips, which increases the stress on 
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Table 1 List of mRNAs with increased polysome association during ER stress (treatment with DTT or TM) 
DTT TM mRNA 

Ribosomal proteins/translation 



2.81 


2.92 


60S ribosomal protein L27a (AFUA_3G05600) 


2.15 


3.20 


40S ribosomal protein S29 (AFUA_6G12720) 


2.87 


2.26 


40S ribosomal protein S8 (AFUA_6G07360)" 


2.80 


2.16 


50S ribosomal protein L36 (AFUA_4G12810) 


2.59 


2.32 


37S ribosomal protein S16 (AFUA_5G08350) 


1.48 


2.60 


Mitochondrial ribosomal protein L1 1 (AFUA_5G1 1830) 


2.01 


1.93 


Ribosome biogenesis protein (AFL)A_8G04790) 


1.18 


2.49 


37S ribosomal protein S5 (AFUA_5G1 1540) 


1.51 


1.93 


40S ribosomal protein S9 (AFUA_3G06970) # 


1.01 


2.42 


40S ribosomal protein Rps16 (AFUA_2G10500)* 


1.07 


2.20 


60S ribosomal protein L6 (AFUA_6G09060) 


4.51 


2.08 


Aconitate hydratase/Mitochondrial ribosomal protein subunit L49 (AFUA_3G08080)' 


3.11 


1.36 


Eukaryotic translation initiation factor 3 subunit elF-Ca (AFUA_1G05200) 


1.72 


2.49 


Translation elongation factor eEF-1B gamma subunit, putative (AFUA_1G17120) 


1.27 


2.52 


GTP binding protein Gufl (AFUA_3G 14350) 


1.86 


1.05 


Eukaryotic translation initiation factor 3 subunit elF-Ck (AFUA_3G09280) 


Cell membrane/cell wall 






1.08 


6.34 


Squalene monooxygenase Erg 1 (AFUA_5G07780) 


1.61 


1.43 


Ergosterol biosynthesis protein Erg28 (AFUA_2G1 1550) 


1.37 


1.39 


1,3-beta-glucanosyltransferase Gel2 (AFUA_6G1 1390) 


3.51 


4.88 


Cell wall proline rich protein (AFUAJG 13450) 


2.59 


3.62 


Actin cortical patch protein Sur7 (AFUA_2G02310) 


4.21 


2.23 


CFEM domain protein (AFUA_6G 14090) 


1.36 


1.89 


GPI anchored protein (AFUA_2G07800) 


1.78 


1.01 


GPI anchored dioxygenase (AFUA_3G01 800) 


1.24 


3.13 


Dolichol phosphate-mannose biosynthesis regulatory protein Dpm2 (AFUA_1G03K 


1.32 


1.05 


N-acetylglucosaminyl-phosphatidylinositol deacetylase, putative (AFUA_5G12550) 


1.10 


2.58 


Integral membrane protein (Pthll) (AFUA_6G03600) # 


Protein folding & modification 






1.29 


2.48 


Alpha-1,2-mannosyltransferase (Alg2) (AFUA_5G13210) 


1.83 


2.58 


Disulfide isomerase (TigA) (AFUA_5G 12260)" 


1.04 


2.84 


Protein disulfide isomerase Pdil (AFUA_2G06150) 


1.04 


1.82 


N-acetyltransferase family protein (AFUA_4G10930) 


1.35 


1.47 


N-acetyltransferase complex ARD1 subunit (AFUAJG09600) 


2.11 


2.89 


Prefoldin subunit 5 (AFUAJG 10740)" 


Endosome/protein transport and sorting 






1.71 


1.94 


Rho GTPase activator (Bem3) (AFUA_6G06400) 


1.69 


1.41 


Fasciclin domain family protein (AFUA_1G 14300) 


3.75 


1.50 


Ras-like GTP-binding protein (AFUA_4G03100) 


1.39 


2.26 


Endosomal cargo receptor (Erv14) (AFUA_6G07290) 


1.58 


2.16 


Synaptobrevin-like protein Sybil (AFUA_6G1 1270) 


1.57 


3.04 


RAB GTPase Vps21/Ypt51 (AFUA_3G 10740)" 


5.38 


2.13 


Vacuolar protein sorting 55 superfamily (AFUA_6G04780)" 
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Table 1 List of mRNAs with increased polysome association during ER stress (treatment with DTT or TM) (Continued) 



1.03 


2.77 


AP-1 adaptor complex subunit Sigma (AFUA_2G01570) 


1.93 


1.75 


Mitochondrial import inner membrane translocase subunit (TIM22) (AFUA_5G02200)* 


Transcription 






3.07 


3.04 


bZIP transcription factor JlbA/IDI-4 (AFUA_5G01650)* 


1.88 


3.90 


CBF/NF-Y family transcription factor (AFUA_2G14250) 


4.34 


1.41 


C6 transcription factor (AFUA_3G09130) 


2.03 


3.35 


Transcription factor RfeF (AFUA_4G10200) 


1.03 


4.24 


CP2 transcription factor (AFUAJG17350) 


3.09 


1.71 


bZIP transcription factor (LziP) (AFUAJG16460)* 


1.84 


2.39 


Transcription initiation factor TFIID, 31 kd subunit, putative (AFUAJG14600) 


T 1 Q 

z. 1 o 




rrariscri piion idcior iAruM_oui izjuj 


1.48 


1.08 


RNA polymerase II mediator complex component Srb8, putative (AFUA_3G06250) 


1.63 


3.90 


CHCH domain protein (AFUA_3G06370)* 


2.44 


1.48 


Nitrogen metabolite repression regulator NmrA (AFUA_5G02920) 


Stress response 






1.24 


1.99 


General stress response phosphoprotein phosphatase Psr1/2 (AFUA_1G04790) # 


4.37 


3.43 


Glutathione peroxidase Hyr1 (AFUA_3G12270) 


Values represent log2[translational efficiency ratio], 


as described 


in Methods. # Those mRNAs that were also subject to translational upregulation in the thermal 



stress dataset (37°C, 60 min). 



the secretory pathway [7]. In this study, we found that 
the translational efficiency ratios for ergl and erg28 
mRNAs, encoding two key steps in the ergosterol path- 
way, increased in the presence of DTT or TM (Table 1). 
We have previously shown that these mRNAs do not 
increase in unfractionated total RNA under these same 
conditions [6], suggesting that they are under trans- 
lational control. Northern blot analysis of polysome- 
fractionated RNA confirmed these findings; the amount 
of ergl mRNA that was loaded by 5 or more ribosomes 
was greater in the presence of DTT, indicating enhanced 
translation efficiency (Figure 3). It interesting to note 
that in addition to its enzymatic function, Erg28 has a 
unique role as an ER transmembrane protein that acts 
as a scaffold to tether other members of the ergosterol 
biosynthetic complex into a single functioning unit [34]. 
Thus, an increase in the translational efficiency of erg28, 
and potentially other ergosterol biosynthetic mRNAs, 
could work in concert with UPR-mediated transcrip- 
tional increases to drive flux through the sterol pathway 
and support membrane homeostasis. To our knowledge, 
this is the first evidence that mRNAs encoding ergos- 
terol biosynthetic enzymes are subject to translational 
control in A. fumigatus. Since overexpression of mRNAs 
involved in sterol biosynthesis is an established mecha- 
nism of triazole antifungal drug resistance [35], it is 
intriguing to speculate that an increase in the transla- 
tional efficiency of a mRNA in this pathway, even with- 
out a change in mRNA abundance, could represent 



a previously overlooked mechanism of antifungal drug 
resistance. 

A. fumigatus [3(l-3)glucanoxyltransferases (Gell and 
Gel2) catalyze the elongation of (3(1-3) glucan side chains 
and influence morphogenesis and virulence [36,37]. A pre- 
vious report indicates that both Gell and Gel2 are consti- 
tutively transcribed in A. fumigatus [37]. However, here 
we demonstrate that the translational efficiency of the gel2 
mRNA increases 2.5 fold during ER stress, suggesting that 
an increase in Gel2 protein is needed to protect the wall 
under these conditions. Gel2 contains a glycosylphospha- 
tidylinositol (GPI) anchor that tethers it to the plasma 
membrane [37], which facilitates its role in maintaining 
cell wall integrity. Interestingly, at least three other 
mRNAs encoding GPI-anchored proteins of unknown 
function also showed increased ribosome occupancy dur- 
ing ER stress. In addition, ER stress caused increased poly- 
some association of the mRNA encoding the major 
regulatory component for the rate-limiting step in GPI 
anchor biosynthesis, Dpm2, as well as the subsequent 
enzyme in the pathway, AfPIG-L. Together, these findings 
argue that rapid translation of GPI-anchored proteins is 
necessary to protect the fungus under conditions that dis- 
rupt ER homeostasis, mostly likely due to their role in 
maintaining the cell wall [37-39]. It is worth noting that 
GPI anchor biosynthesis is an emerging target for the 
development of new antifungal therapy [40-42]. Further 
understanding of the mechanism(s) by which translational 
regulation impacts GPI anchor production could suggest 
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Figure 3 The erg 7 mRNA increases its association with 
polysomes during ER stress. Mycelial extracts from control 
(untreated) and TM-treated cultures were fractionated into 7 pools. The 
RNA in each pool was then separated by RNA gel electrophoresis and 
the level of ergt mRNA in each fraction was determined by hybridization 
to an erg1 probe. Band intensities were quantified by phosphorimager 
analysis and shown on the top graph. A representative OD 254 profile is 
superimposed on the graph for reference. The findings demonstrate 
increased erg? mRNA levels in the polysome fraction during ER stress. 



novel strategies to enhance pharmacologic inhibition of 
this pathway. 

Host-temperature adaptation involves distinct 
translatome remodeling 

The primary ecological niche for A. fumigatus in nature 
is composting organic material, an environment that 
undergoes constant fluctuations in temperature as a 
consequence of complex microbial activity. A. fumigatus 
has evolved mechanisms to thrive under these condi- 
tions, and is considered one of the most thermotolerant 
species of mold [43]. Since elevated temperatures induce 
conformational changes in proteins [44], an increase in 
temperature is likely to engage pathways that are rele- 
vant to ER stress response. We therefore compared the 
translational efficiency of A. fumigatus mRNAs at 25°C, 
representing the environment, to that of mRNAs follow- 
ing a shift to 37°C, reflecting adaptation to the mam- 
malian host. Ribosome fractionation showed that total 



polysome levels increased within 30 min of the shift to 
37°C, consistent with the need for increased proteins at 
this optimal growth temperature (Figure 4). Polysome 
peak heights declined somewhat after 60 min at 37°C, 
presumably reflecting a return to steady-state levels at 
the new temperature. Two criteria were employed to 
define differentially translated mRNAs during this tran- 
sition. First, we considered all mRNAs that shifted from 
fraction-U to fraction-W following the temperature shift 
to have a temperature-induced increase in translational 
efficiency (two-fold cutoff). This resulted in the identifi- 
cation of 311 translationally upregulated mRNAs 30 min 
after the temperature shift, and a total of 499 mRNAs at 
the 1 h time-point. Some of these mRNAs could also be 
upregulated at the level of transcript abundance during 
ER stress. Thus, in order to enrich for mRNAs that are 
predominantly regulated at the level of translational 
efficiency, the dataset was narrowed to those mRNAs 
that showed a minimal two-fold increase in translational 
efficiency ratio when normalized to relative transcript 
abundance in unfractionated RNA. Applying these cri- 
teria, 78% of mRNAs were translationally upregulated at 
the 30 min time-point and 75% were upregulated at the 
1 h time-point. These findings demonstrate that thermal 
stress is similar to DTT- and TM-induced ER stress 
in its reliance on translational regulation as a rapid- 
response mechanism to manufacture essential proteins 
that are needed to protect the fungus during host- 
temperature adaptation. 

Hierarchical clustering of all mRNAs that showed 
temperature-dependent increases in translational efficiency 
fell into three major clusters (Figure 5). The first group 
('early') showed a transient increase in translational effi- 
ciency at 30 min that returned to baseline levels by 1 h. 
The second group ('late') showed baseline levels at 30 min 
but an increase at 1 h. The third group ('continuous') 
showed an increase at 30 min that was sustained at 1 h or 
subject to a further increase. Over- represented functional 
groups in the entire dataset of translationally upregulated 
mRNAs at 37°C included nucleotide metabolism (28), 
ribosome function (18), oxidative phosphorylation (26), 
TCA cycle (8), cell cycle (23), and secondary metabolism 
(18) (Additional file 3). The increased translation of 
mRNAs encoding proteins with roles in metabolism fol- 
lowing the temperature shift is consistent with the fact that 
A. fumigatus grows more rapidly at 37°C than it does at 
25°C. However, some metabolic genes were also enriched 
in the downregulated category (see the full dataset, 
ArrayExpress accession E-MTAB-2027), indicating that 
complex metabolic adjustments are operational during the 
transition from 25°C to 37°C. Interestingly, we found that 
mRNAs encoding heat-shock proteins were largely absent 
from the dataset of translationally upregulated mRNAs 
following the shift from 30°C to 37°C. However, this is 
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Figure 4 Effects of a temperature shift on the polysome profile A. fumigatus. Overnight cultures at 25°C were shifted to 37°C and polysome 
profiles were prepared from mycelial extracts at the indicated times. 



congruent with previous transcriptomic analyses, which 
showed little upregulation of heat-shock protein-encoding 
mRNAs under these conditions [45,46]. We speculate 
that the transition from 30°C to 37°C is a relatively mild 
stress for this thermotolerant fungus, requiring minimal 



upregulation of heat-shock proteins for protection, either 
at the transcriptional or translational levels. 

The Venn diagram in Figure 6 illustrates the overlap 
between the transcripts that were translationally upregu- 
lated in response to the 3 types of ER stress applied in 
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Figure 5 Translational response to thermal stress. Hierarchical clustering of mRNAs that showed temperature-dependent increases in translational 
efficiency fell into three major clusters (left). A schematic representation of each pattern of translational efficiency change is shown in the center: Late 
(increased at 60 min), Continuous (increased at 30 min and maintained or further increased at 60 min, and Early (transient increase at 30 min). The Venn 
diagram illustrates the number of mRNAs in each category. 
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Thermal shift 

Figure 6 Summary of translationally upregulated mRNAs. The 

Venn diagram shows the number of transcripts with a higher 

translational efficiency (2-fold change) when treated with DTTorTM 

or when moved from 25°C to 37°C for 1 h. The overlapping region 

(shaded) represents a core set of translationally up-regulated genes 

that respond to all three types of stress in A. fumigatus. 
\ ) 

this study. A total of 44 mRNAs were shared between 
the three datasets, only one of which encoded a classical 
ER stress response protein involved in protein folding, 
TigA (Additional file 4). These findings demonstrate 
the remarkable breadth of functions that underpin ER 
homeostatic functions in this organism. 

The ER stress translatome contains a UPR-dependent 
mRNA isoform 

The regulation of mRNA translation efficiency can occur 
by a variety of mechanisms, one of which includes alter- 
ations in the mRNA sequence that enhance translatability 
[47]. To determine whether ER stress induces changes in 
mRNA structure, whole transcriptome shotgun sequen- 
cing (RNA-seq) was performed in the presence or absence 
of DTT. The RNA-seq coverage plot for the hacA mRNA 
served as a control for this analysis, revealing splicing of 
the same conventional and unconventional introns that 
have been validated experimentally [7] (Figure 7). A 
substantial amount of unconventional intron splicing was 
evident in the absence of DTT, consistent with our pre- 
vious finding that A. fumigatus relies upon the UPR even 
during normal vegetative growth [6]. This contrasts the 
situation in S. cerevisiae, where there is minimal splicing 
of HAC1 in the absence of exogenous ER stress [27], pos- 
sibly due to the decreased demands on the secretory path- 
way in a unicellular yeast. As expected, treatment of 
A. fumigatus with DTT was associated with an increase in 
reads across the length of the hacA gene (Figure 7), with 
the exception of the unconventional intron, demonstrating 



an increase in the ratio of spliced to unspliced hacA RNA 
under these conditions. The unconventional intron in 
A. fumigatus hacA is only 20 nucleotides in length, which 
is similar to what has been reported in the human homo- 
log but strikingly different from the 252 nucleotide intron 
in the S. cerevisiae homolog. In S. cerevisiae, the un- 
conventional intron blocks translation of the mRNA by 
forming a stem-loop structure with the 5'UTR [48]. The 
removal of the intron by Irel-mediated splicing releases 
this translation block, allowing the spliced mRNA to be 
translated. The small size of the hacA intron in A. fumiga- 
tus makes a similar translation block mechanism unlikely, 
similar to what has been reported in mammals, Caenor- 
habditis elegans, Candida albicans, and other filamentous 
fungi [12,49-52]. In fact, the unspliced mRNA in humans 
is translated into a protein product that contains a hy- 
drophobic segment that tethers the mRNA to the ER 
membrane, thereby facilitating splicing by Irel [53]. In 
A. fumigatus, both the unspliced and spliced hacA mRNAs 
can be readily identified in fraction-W by RT-PCR (data 
not shown), suggesting the possibility that the unspliced 
RNA is translated. It will be interesting to determine 
whether its putative encoded product is involved in a simi- 
lar ER membrane tethering mechanism in A. fumigatus. 

We next analyzed the RNA-seq profiles of all 233 trans- 
lationally upregulated mRNAs identified in our ER stress 
study (Figure 2). The RNA-seq coverage plot of the 
mRNA encoded by AfuA_3G 13490 showed a striking 
change in the presence of DTT (Figure 7). This mRNA 
encodes the A. fumigatus homolog of yeast Yvcl, a transi- 
ent receptor potential (TRP) channel protein in the va- 
cuolar membrane that is the major release mechanism for 
intracellular calcium stores [54]. In the absence of DTT, 
the number of sequence reads was comparable along the 
length of the yvcl mRNA (Figure 7, red tracing), with the 
exception of 4 predicted introns denoted by the vertical 
columns. However, DTT treatment induced an increase 
in sequence reads, but only at the 3'-end of the gene 
(Figure 7, blue tracing). This mRNA did not splice out in- 
trons 3 and 4, suggesting that DTT stress was inducing a 
novel mRNA isoform derived from the yvcl transcription 
unit, henceforth referred to as yvcla. Northern blot ana- 
lysis using the full-length yvcl open reading frame (orf ) as 
a probe confirmed that ER stress induced yvcla expres- 
sion, but osmotic stress with NaCl did not (Figure 8). In 
addition, DTT failed to induce yvcla in two UPR mutants, 
AireA and AhacA, indicating that its presence is both ER 
stress-specific and downstream of the UPR. Sequence ana- 
lysis of the yvcla cDNA identified a single long open 
reading frame that would encode the C-terminal 127 
amino acids of the full-length Yvcl protein (accession #: 
XP_001481630.1). Although the oligonucleotide used for 
microarray hybridization would not distinguish yvcla 
from yvcl, RT-PCR analysis confirmed that both mRNAs 
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AfuA_3G04070, hacA AfuA_3G13490, yvcl 

Figure 7 RNA-seq coverage plots for the hacA and yvd mRNAs. The number of sequence reads on the y-axis (reads per kilobase per million) 
is shown along the length of each gene in the absence (red) or presence (blue) of ER stress (1 mM DTT, 1 h). Vertical lines demarcate predicted 
intron boundaries (shown in green for the unconventional intron in hacA). The coverage plot for yvcl shows an increase in reads at the 3' end of 
the gene specifically in the presence of ER stress. 



are located in fraction-W during ER stress (data not 
shown), suggesting that both of them contribute to the ER 
stress response. No alternative YVC1 isoforms could be 
identified in S. cerevisiae treated with either DTT or TM 
(data not shown), suggesting a fundamental difference be- 
tween fungal genera that may reflect the increased com- 
plexity of ER stress mechanisms among filamentous fungi. 
Since yeast Yvcl has a role in calcium homeostasis [54], it 
is appealing to speculate that the alternative isoform of 
this protein provides a novel regulatory mechanism to link 
ER stress with calcium signaling. The precise nature of the 
protein product encoded by this novel mRNA, its contri- 
bution to the biology of A. fumigatus, and the molecular 
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Figure 8 Northern blot analysis of yvcl expression. RNA was 

isolated from the indicated strains in the presence or absence of ER 
stress (1 mM DTT for 1 h, 10 ug/ml TM for 1 h) or osmotic stress 
(0.8 M NaCI for 1 h). The blot was hybridized with the full-length 
A. fumigatus yvcl gene, then stripped and re-probed with the UPR 
marker bipA. Accuracy of RNA loading is shown by SYBR-green 
staining of rRNA. 



mechanisms that enhance its ribosome occupancy during 
ER stress are currently under active investigation. 

Conclusions 

In summary, this study reports the first global analysis 
of transcript-specific translational regulation during ER 
stress in the pathogenic mold A. fumigatus. The results 
define a core ER stress translatome and demonstrate 
that translational regulation is a major component of the 
response of A. fumigatus to environmental conditions 
that perturb ER homeostasis. We also provide evidence 
that the ER stress translatome contains a previously 
unidentified mRNA that is induced by ER stress in a 
UPR-dependent manner. Together with our previous 
analysis of the UPR transcriptome, these findings begin 
to develop a comprehensive understanding of how 
A. fumigatus responds to ER stress. Since ER stress is 
experienced by several fungal pathogens in the host 
environment [6-10], further understanding of these 
pathways, and the mechanisms used to deploy them, 
may provide a new perspective on fungal pathogenesis 
and offer novel strategies for therapeutic intervention. 

Methods 

Strains and culture conditions 

The wt strain used in this study is AfS28 (AakuAwptrA) 
[55]. Conidia were harvested from colonies grown on 
OSM plates {Aspergillus minimal medium containing 
10 mM ammonium tartrate and osmotically stabilized with 
1.2 M sorbitol). For analysis of the ER stress response, a 
250 ml flask containing 50 ml of YG medium (0.5% yeast 
extract, 2% glucose) was inoculated with 5 X 10 7 conidia 
and incubated for 16 h with shaking (200 rpm) at 37°C. ER 
stress was then induced by treating the cultures with 
1 mM DTT or 10 ug/ml TM for 1 h. Hyperosmotic stress 
was induced by adding NaCI to a final concentration of 
0.8 M and continuing the incubation for an additional 
hour. For analysis of the heat-shock response, 250 ml 
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flasks containing 50 ml of YG medium were inoculated 
with 5 x 10 conidia and incubated for 16 h with shaking 
(200 rpm) at room temperature. Thermal stress was then 
applied by transferring the flasks to a shaking 37°C incuba- 
tor (200 rpm) and harvesting the mycelium after 30 min 
and 60 min of incubation. 

Polysome fractionation and RNA extraction 

Sample preparation and polysome analysis was performed 
as previously described, with modifications [56] . Following 
the induction of ER or thermal stress in the cultures, 
translating polyribosomes were halted on the mRNAs by 
the addition of cycloheximide to 0.1 mg/ml and incuba- 
ting at 37°C for an additional 5 min. The culture was then 
chilled in an ice bath for 5 min prior to harvesting the 
mycelium. The hyphae were washed twice with 5 ml of 
lysis buffer (10 mM Tris-HCl [pH 7.5], 100 mM NaCl, 
30 mM MgCl 2 , 0.1 mg/ml cycloheximide, and 0.2 mg/ml 
heparin), flash frozen in liquid nitrogen and mechanically 
crushed. After resuspending in 0.5 ml of lysis buffer, the 
lysate was cleared by two subsequent microcentrifugation 
steps (15,000 x g, for 5 min at 4°C) and the RNA content 
in the supernatant was quantified by absorbance at 
260 nm. Equal amounts of RNA (20-30 A 26 o units) were 
loaded onto a 12-ml linear sucrose gradient (7 - 47%) pre- 
pared in gradient buffer (50 mM Tris-acetate, 50 mM 
NH 4 C1, 12 mM MgCl 2 , 1 mM DTT, and 0.2 mg/ml hep- 
arin). The gradients were centrifuged at 150,000 x g for 
2.5 h at 4°C, using a Sorvall SW 41Ti rotor. Gradient ana- 
lysis was performed using an ISCO gradient collector with 
continuous monitoring at 254 nm. Individual fractions 
were collected with a Foxy Jr. fraction collector and RNA 
was precipitated from 0.5 ml fractions by mixing with an 
equal volume of 6 M guanidine thiocyanate and 2 volumes 
of 100% ethanol and incubating overnight at -20°C. The 
RNA was pelleted, washed and resuspended using stand- 
ard procedures. For microarray analysis, RNA from frac- 
tions containing less than 5 ribosomes/ mRNA ('U') or 5 
or more ribosomes/mRNA (*W) were pooled and pre- 
cipitated with 1.5 M LiCl, followed by washing to remove 
residual heparin. For northern blot analysis of ergl expres- 
sion, the sucrose gradient was divided into seven sequen- 
tial fractions representing the entire gradient, and the 
RNA was precipitated as indicated above. For experiments 
that required unfractionated RNA (unfractionated 
controls for the thermal shift microarray experiment, 
northern blot analysis of ergl mRNA, and RNA-seq ana- 
lysis of DTT-treated cultures), the mycelium was crushed 
in liquid nitrogen and total RNA was extracted using the 
TRIZOL method [57]. 

Microarray hybridization 

The RNA labeling reactions and hybridizations were 
performed as described in the J. Craig Venter Institute 



(JCVI) standard operating procedure http://pfgrc.jcvi. 
org/index.php/microarray/protocols.html) and transcrip- 
tional profiles were generated by interrogating the Af293 
spotted oligonucleotide microarray containing 10, 503 
spots. Each gene was present in triplicate on the array, 
and all hybridizations were repeated in dye swap experi- 
ments. The data for each gene were averaged from the 
triplicate genes on each array and the duplicate dye swap 
experiment (a total of six readings for each gene) and 
the gene expression ratios were log2-transformed. Plot- 
ting open reading frame length against fold increase in 
the W fraction showed no bias towards longer tran- 
scripts, indicating that an increase in ribosome loading 
on a particular transcript is not an artifact of mRNA 
length (data not shown). 

Functional annotation of genes present within the dataset 
was analyzed using FungiFun [58] and enrichment of func- 
tional groups was performed using FunCat method. Hier- 
archical clustering was performed using Cluster 3.0 [59] 
and the cluster tree was visualized using JAVA Treeview 
[60]. All RNA samples were hybridized with a reference 
sample obtained from Af293 in order to allow for cross- 
comparison. The translational efficiency of individual 
mRNAs during DTT/TM treatment was defined as the 
ratio of the hybridization signal in fraction-W over that of 
fraction-U, using a 2-fold difference between conditions as 
the cut-off value for a change in translational efficiency. 
Normalization to total mRNA abundance was not per- 
formed because the mRNAs that fit these criteria showed 
no increase in abundance under the same conditions [6]. 

The translational efficiency of individual mRNAs at 25°C 
and following a temperature shift to 37°C (after 30 min or 
60 min) was defined as the ratio of the hybridization signal 
in fraction-W over that of fraction-U, using a 2-fold 
change between conditions as the cut-off value for a 
change in translational efficiency. In order to enrich for 
mRNAs that are predominantly regulated by changes in 
translational efficiency (as opposed to transcript abun- 
dance), the dataset was normalized to transcript levels in 
unfractionated RNA. RNA abundance was determined by 
interrogating the microarrays with unfractionated RNA 
and the change in the translational efficiency of each 
mRNA upon thermal shift was calculated as (fraction-W/ 
fraction-U) /total transcript abundance. 

RNA sequencing 

RNA-seq was performed by the Genomics Sequencing 
Core (GSC) at the University of Cincinnati. Using TruSeq 
RNA sample preparation kit (Illumina), total RNA (RIN > 
7.0, Agilent 2100 Bioanalyzer) was converted into a library 
of template molecules suitable for subsequent cluster 
generation and sequencing by Illumina HiSeq. Poly(A) n 
mRNA was extracted and fragmented into smaller pieces 
(-140 nt). The cleaved RNA fragments were converted 
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into first strand cDNA using reverse transcriptase and ran- 
dom primers, followed by second strand synthesis using 
DNA polymerase I and RNAse H. The cDNA fragments 
were then subject to end-repair followed by addition of a 
single A base and ligation of adapters. The products were 
indexed individually, purified and enriched by PCR to 
create the final cDNA library. The generated library was 
validated and quantified using Kapa Library Quantification 
kit (Kapabiosystem). Six individually indexed cDNA libra- 
ries of equal amounts were pooled for clustering in cBot 
system (Illumina). Libraries were clustered onto a flow cell 
using Alumina's TruSeq SR Cluster Kit v3, and sequenced 
for 50 cycles using TruSeq SBS kit on Illumina HiSeq 
system. 

FASTQ files containing 50 bp single-end RNA-Seq 
reads were mapped to the Aspergillus fumigatus genome 
sequence (taxid:330879) by TopHat [61]. Transcript as- 
sembly and abundance estimation were performed by 
Cufflinks [62]. 

Reads corresponding to 233 genes of interest were 
filtered and the coverage of each nucleotide position was 
counted using a semi-automated method in order to 
ensure accuracy of analysis. Coverage plots for each of 
the 233 genes under two conditions were plotted using 
Matlab®. 

Analysis of mRNA expression by northern blot analysis 
and qPCR 

RNA samples were fractionated by formaldehyde gel elec- 
trophoresis, and visualized by SYBR green staining. The 
RNA was then transferred to BioBond nylon membranes 
(Sigma) and hybridized to a 32 P-labeled DNA probe as 
previously described [7]. Probes specific for the A. fumiga- 
tus ergl, yvcl and bipA genes were PCR amplified from 
genomic DNA using the following primers: ergl: 5 - 
CGTCAGTGTTGTTGAGAC-3 ' and 5'- GAAGGTCGA 
GAGCTGCTTC-3 '; yvcl: 5'- CAATGCTGTGGACGA 
GTACATG-3 ' and 5' - GTGCTCCTCTGTATCCTTC 
TTC-3'; bipA: 5'- GTCTGATTGGACGCAAGTTC-3 ' 
and 5 - ATCTGGGAAGACAGAGTACG-3 ' . Hybridi- 
zation intensities were quantified by phosphorimager ana- 
lysis using Image Lab software. 

For qPCR analysis one ug of RNA from pooled fractions 
corresponding to fraction-U or fraction-W was reverse- 
transcribed with M-MuLV reverse transcriptase (NEB) 
using oligo (dT)i 8 and 18S rRNA primers (primer 
713-TGAGCCGATAGTCCCCCTAA and primer 714- 
GACTCAACACGGGGAAACTC). The qPCR was per- 
formed using the iTaq™ universal SYBR* green supermix 
(Bio-Rad) according to the manufacturer's protocol. The 
melting curve was monitored to verify specificity of the 
amplification reaction. Controls reactions in the absence 
of reverse transcriptase were used verify the absence 
of DNA contamination. The 18S rRNA present within 



fraction-U or fraction-W was used as an endogenous con- 
trol to derive a ACt value for each fraction. A translational 
efficiency ratio (W/U) was derived by subtracting ACt of 
fraction-W from that of fraction-U, representing AACt. 
Change in W/U ratios upon treatment with DTT or TM 
was then plotted using 2-AACt of untreated samples as 
the reference. Primers used for qRT-PCR are as follows: 
(3-tubulin (AfuA_lgl0910), primer 554-CACGGATCTT 
GGAGATC and primer 562-ACAACTTCGTCTTCGG 
CCAG; squalene monooxygenase ergl (AfuA_5g07780), 
primer 810-AGCTGCGATCTATGCCGAATTCCT and 
primer 799-TCCCAGTTGGAAGTAACGGAAGCA; vacu- 
olar protein sorting 55 superfamily vpsSS (AFUA 6G04 
780), primer 804-GCGCTCTCCTTTGTTCTTGCCATT 
and primer 805-AAGACCTCCGAGGATGGACATGAT; 
bZIP transcription factor jlbA/IDI-4 (AFUA_5G01650), 
primer 813-TTGATGTGAACGACTCTCTGCCGT and 
primer 814-TAGCTTCGACACCCGCATCTTCAA. The 
data were compared by Student's t-test and a p < 0.05 was 
considered significant (indicated by the asterisk, Additional 
file 1). 

Data availability 

The microarray and RNA-seq data sets reported in 
this article are available in the ArrayExpress database 
(microarray accession E-MTAB-2027, RNA-seq acces- 
sion ERP004296). 

Additional files 



Additional file 1: Validation of the translationally regulated dataset 
by qPCR. The levels of 1 8S rRNA in fraction-U or fraction-W was used as an 
endogenous control to derive a ACt value for each fraction. A translational 
efficiency ratio (W/U) was then calculated by subtracting ACt of fraction-W 
from that of fraction-U, representing AACt. The change in W/U ratios upon 
treatment with DTT or TM was then plotted using 2-AACt of untreated 
samples (UT) as the reference. 

Additional file 2: List of mRNAs with decreased polysome 
association during ER stress (treatment with DTT or TM). Values 
represent log2 [translational state efficiency], as described in Methods. 

Additional file 3: List of over-represented KEGG pathways in the 
dataset of translationally regulated mRNA following a shift to 37°C. 

Additional file 4: List of mRNAs with increased polysome 
association during each of the three forms of ER stress: treatment 
with DTT, TM and thermal stress. Values represent log2[translational 
efficiency ratio], as described in Methods. "mRNAs subject to translational 
upregulation in the thermal stress dataset at 60 min. 
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